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The transition of the laminar 


boundary layer into the turbulent 
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state, which results in an increased drag, can be avoided by sucking 
off the boundary la-ye.r- .par.ti.cLes .near the:- wall..- -,The technically-in- 
teresting case of sucking the particles using individual slits is in- 
vestigated for bodies of revolution in incompressible flow. 

The boundary layer calculation is done using the integral con- 
ditions for momentum and energy, in order to determine the laminar- 
turbulent transition point, a new semi-empirical criterion is intro- 
•duced. The changed boundary layer t quantities behind a suction slit 
are approximately determined using] the remaining profile after the 
suction process. In order to evaluate the slit suction, we carry out 
a drag calculation with consideration-of the suction power. 

The calculation method uses a j ALGOL program. The reliability of 
the method is confirmed by comparison calculations.. 

The results of the variational calculations showt that there is 
an optimum suction height, where the slot separations are maximum. 
Combined with favorable shaping of] the body, it is possible to keep 
the boundary layer over bodies of revolution laminar at high Reynolds 
numb.ers using relatively few suction' slits and small amounts of suc- 
tion flow. 
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1. Notation 


b. 


V > 


8\J } b ^5 
A^, Ajj 


D 


"Q 

d 

F 

'^13 ^ 2 ^ ^4 

.^12" ^32 
•1 

•L 

m 


m- 


■E 


n, N 

n. 


Q 

Q 


-N, 
0 
P 

q 

Q 

'r 

Re 

■s 

t 


coefficients 

coefficients 

amplification ratio 
slot length 
exponent 

dissipation coefficient 

friction coefficient 

pressure coefficient 

suction amount coefficient 

drag coefficient 

equivalent suction drag coef- 
ficient 

maximum thickness 
area 

universal functions 
shape parameter 
length 

Mangier-constant 

Hartree parameter 

Hartree parameter of plane 
replacement flow 

boundary layer variables 

number of suction points 

suction -power 

surface area 

static pressure 

stagnation pressure 

suction yield /6 

suction amount 

body radius 

Reynolds number 

slot v^idth 

time 
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u , u 

<3x qr 

qt 

U, V 

V 

w 

w 

Q 

X, r, 

X’ 

< z 

a, 3 

IT 
X 


9 

z 


6 

✓v 

6. 


'6- 


‘Q 

'IQ 

’2Q 

'3Q 


e 

'n 


^Q* 

j? 

0 

K A 

u 

V 
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9o 


Induced velocity components 

induced tangential velocity 

velocity components 

average velocity in suction slot 

volume 

drag 

equivalent suction drag 
cylinder coordinates 
cartesian coordinates 
coordinate along body contour 
suction height 
' thickness parameter 

boundary layer variables 
buildup variable 

f ' 

cone angle, wedge angle 

I ■ wedge angle of plane replace- 
ment flow 

, boundary layer thickness 

1 

, displacement thickness 

! momentum loss thickness 

; energy loss thickness 

suction loss thickness 

' suction displacement thickness 

' suction momentum loss thickness 

' suction energy loss thickness 

boundary layer variable 

relative wall separation, 

; radius coordinate /7 

, relative suction height 

' ‘^all point angle 

I momentum loss area 

I Pohlhausen-shape parameter 

I 

j dynamic viscosity 

j kinematic viscosity 

i length coordinate 

i density 

I nose radius 
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I 

II 
* 


apex radius 
shear stress. 

I 

i 

1 

( 

■ referred to separation point 

; referred to suction region 

I nose part 

j end of rotationally-symmetric 

; boundary layer calculation 

■ referred to the area 

referred to the velocity trial 
solution 

1 tail part 

i 

I referred to tail part 

' referred to the instability 

I point 

i referred to length 

I ■ laminar 

I middle part 

• referred to the surface 

.referred! to potential flow 
' referred to suction point 

referred to support point 
turbulent 

referred to transition point 

referred to the volume 

referred to the wall 

[ referred to edge of boundary 

layer 

referred to momentum loss 

t ‘ thickness 

( 

' referred! to incident flow 

! ahead of suction slit 

I 

] behind suction slit 

average value 
transformed variable 
! nondimensional quantity. 
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' 2. Introduction . /9 

. Bodies of !revo|lution are primarily used in aerodynamics and hydro- 

> 

•dynamics in order to support useful loads (for example, aircraft bodies., 
underwater vehicles). Therefore, optimum bodies of revolution from the 
I flow point of view are those whosejflow resistance with respect to the 
•useful volume is a minimum. Minimum specific drag is achieved by the 
'following: | 

- by means of a large volume related to the surface in the flow, 

• - by a small drag coefficient [referred to the surface area. 

Whereas the volume-to-surface area [is increased with increasing thick- 
ness ratio, the drag-<coef fi-cient -re-ferred— to--the surface can be reduced 
•'by keeping the boundary layer laminar .in the_ flow. 

In contrast to conventional shapes with a long cylindrical central 
[part and a‘ small thickness ratio, Hertel’Cl to 3] suggested spihdle- 
i shapes for aircraft bodies, which are characterized by a large thick- 
'ness ratio and a pointed nose part^ where the laminar boundary layer is 

' ' i 8 

I stabilized. At high Reynolds numbers (Re^ > 10 ) ^ such as occur in 
J aviation and underwater technology j the laminar effect which can be 
j reached by shaping alone is not Important. Already by sucking off 
, small amounts of flow, the laminar [boundary layer can be additionally 
[stabilized, so that the boundary layer transition .into the turbulent 
t state is delayed, and the drag 'is reduced with considerationHio f e; 
i suction power. The amounts sucked off required to Influence the boun- 
Idary layer are smallest for a continuous distribution over the surface. 

i 

;This optimum case is the simplest to analyze theoretically, but cannot 
, 1 

,be realized in practice ^because j o f the technical .effort for suction. 

! In technical applications, one can only consider sucking off the 

, f low through slits arranged perpendicular to the flow direction, or 
'using gaps or perforated strips. The effectiveness of such suction slits 
I and slots was demonstrated' by .Pfenninger [4] to E6] and Lachmann [73 
Ifor airfoils. Based oh thej large number of suction points, the men- /lO 
tioned suction configurations are not completly satisfactory for tech- 
nical applications. 

In the case of spindle-shaped 1 bodies of revolution, it seems Ipr o-i 
rising to expand the laminar flow range even at high Reynolds numbers, 
using a relatively low number of suction points, so that the drag re- 
' ductions considering the suction power will be worthwhile. ; 
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3. Problem Formulation and Assumptions 

I 

f 

In the present paper^ we will | theoretically investigate the tech- 
nically interesting case of maintaining a boundary layer laminar using 
.suction slits on bodies of revolution. In particular, as a function 
of the body shape and the Incident {Reynolds number, we will determine 
the following: 

- optimum position of suction | slits for minimum number of slits, 

- the flow drag with consideration of suction power. 

The investigations are restricted to the following: 

- bodies of revolution with flow in the axial direction, 

- incompressible flow ; 

- flow without' any disturbances, 

- hydraulically-smooth surface 

The calculation method to be designed includes the following fac- 

1 

tors : 

- calculation of the laminar and turbulent boundary -layer over 

bodies of revolution. with consideration of Influencing the laminar 
(boundary layer. | 

- determination of laminar-turbulent boundary layer transition 
for bodies of revolution. 

- influencing of. laminar bound.ary layer using suction slits. /II 

As will be shown later on, we can neglect the following: 

- pressure losses by suction slits, throttling points and lines 
when estimating the suction power, 

- the sink effect of the suction slits. 

The computer program to be es-t'ablished will allow the following 
systematic variations: 

- variations in body shape i 

- variations in suction condit;ions 

- variations in Incident Reynolds number. 

In addition, the program is to include the case without suction 
■and the calculation of flat contours. 

Based on available measurement* results and theoretical -solutions , 
we will demonstrate the reliability of the calculation 'method. In ad- 

j 

dition, we will carry out several characteristic variational calcula- 
tions . ! 
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■ 4, Contour and Pressure Distribution of Bodies of Revolution 
,4.1, Description 'of Body Contour 


I 

One can only vary the body shape in a systematic way if the body 
lOutline can be described mathematically. 

( A -body of revolution consists {of a nose part, a central part, and 

!a tall part, as- shown i-n Figure 1. ' In most cases the central part is 

* j 

cylindrical. In the ease of the spindle shapes described in [2] and- 
,[3], the lengjths' of the central part and the nose radius are zero. 

The mathematical description of the body contours is the simplest /12 
‘if separate trail solutions are assumed for the nose, central, and tall 
•parts. In this ca^s'ej'the 'thlc'Kness' and' lengt'h" ratios do not occur in 

t • j • 

,the contour functions, but only determine their coordinate system. For 
i \ 

tthe other dimensions shown in Figure 1, shaped parameters are introduced 
which are independent of the thickness and length ratios. 

* If the coordinates for the nose, central, and tail- contours are 

‘ t 

; specified as shown -in Figure 1, the;n the following relationships exist 
Sbeteen the body of revolution outlined r(x) and the contour functions • 






i-- 


Nose 


1^/1 - x /1 

"X7T- 


r /I 
B ' 


Central Part 

5m = 


M/I - l^/l 


7T 


1b ■ 


% 


v/' 


■d72T 


5„ 


1 h 


Tail 

'x/r-Ti-r/iF 

n 


iTTT 




(4.1) 


4 /a 


' The requirements for the contour functions 

» 

1 - consideration of the most important body shape parameters 

I 

- and large shape selection | 

_’-for the- most part satisfy the polynomial trial solutions of Koschmieder 
and Walz [8]. In addition, if ther^e is no central part, there can be 
"a transition between the nose part and the tail part without any jump 
in curvature . 

t 

• For body shapes with inflection points, the basis functions of 
,'Oehler [9] are especially suitable.^ By superimposing these basis func- 
^tions, a large selection of shapes -results. 


1.0 



4 >2. Remarks Regarding the Pressure Dir'stributi'on 

In the case of a body -with boundary layer suction, the pressure 
distribution caused. by the body shape is superimposed with a sink ac- 
tion of the suction slits. 

Except for the immediate vicinity of the tail tip, it is possible /13 
to determine the •p:’essure distribution of a body of revolution caused by 
the shape 

c = 1 


f(^) 


"00 


C4.2) 


using potential thepry... _ ^ -j.. 

In this' investigation, we will not discuss the calculation of the 
potential theory pressure distribution caused by the shape. Instead, 
we will indicate the calculation procedure of Oehler [9]> which is es- 
pecially suited for thick spindle shapes with high underpressure peaks, 
The s in 
Chapter 7-4. 


The sink- action of ring-shaped suction slots will be discussed in 




I 5. Boundary Layer Calculation and Drag Layer Calculation of Bodl'es 
i of Revolution ' 

: ■ I 

5.1 Preliminary Remarks I 

! In the case of calculation of 1 the boundary layer over bodies of 

, ■ I 

! revolution with ring-shaped suction slits, one must disbinguish the 
J following: 

; - boundary layer development along the impermeable wall, 

; - influencing the boundary layer at the suction slits. 

! In Chapter 7 we will discuss influencing the laminar boundary 

layer using suction slits. 

• Schlichting [10] gave a summary about boundary layer theory. The 

Prandtl boundary layer equation is the point of departure for the -boun- 
dary layer calculation. One must distinguish between the following: 

- - the exact solutions which satisfy the bo.undary layer equation 
at any point, and 

'■ - the approximate solutions which on the average satisfy the boun- 

dary layer equation over the boundary layer thickness » 

It is only possible to exactly calculate the boundary layer equa- /l4 

i . * 

tlon- for- laminar .b-oundary J.ay-er.s... „ In.-±his.^ -case.,.. .there _ are._-S.uhsi;antial 
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difficulties for arbitrary body shapes and pressure -distributions . The 
approximation theory can be used successfully in most practical problems , 
as ¥alz [11] showed, especially. 

In the case of rotationally symmetric flow, the calculation of 
the boxmdary layer is more difficult than for .plane flow. This is be-' 

; cause the pressure distribution and the body shape appear in the boun- 
dary layer equation. The relationship between the plane boundary layer 
•and the rotationally symmetric boundary layer was analyzed by Mangier 
•[ 12 ]. 

! The similar solutions of the Hartree [13] boundary layer equation 

dete rminedl by Hartre.e._[,13-]-.ar-e the basis ,for_.dey:e.loRing the approximate 
solutions. The most recent trial solutions for the laminar boundary 

r 

Jvelocity profile of ¥alz, Il4], Magnler [15], Geropp [l6] ,' attempt to 
I give the best possible approximation to the Hartree profile. The trial 
I solutions for the turbulent veloci-fey profile are empirical, 
j The appro|ximate [solution principle consists of satisfying the 

•integral and boundary conditions for a velocity profile which is as- 
I sumed to be known. These are deri-^ed from the boundary layer equation. 

I The accuracy of -the approximate method increases with the number of si— 

! multaneously satisfied . boundary layer equations. 

' A quadrature method can- be used to exactly calculate the laminar 

I boundary layer for a. pressure drop,' which is based on the integral 
■•dition for momentum and the wall adhesion. 


/15 


I As Walz [ 17 ] showed later on, jby using a [single^ parameter method 

’(with a shape parameter of the velocity profile), it is possible to 
* obtain results which are more reliabie,' , if instead of the wall condi- 
*tion, one introduces the energy integral condition. 

In place of turbulent boundary layers, the accuracy which can be 
'achieved depends on the empirical law for wall [shear] stress and dissi- 
pation. A corresponding comparison can be found in Pernholz [l8], 

» Just behind a suction point there are laminar velocity points, 

.which differ substantially from the Hartree profiles. The authors 
‘Igli-sch. £ 19 ] and Schllchting, Bussiiiann [20] found similar solutions for 
the suction profiles df the plane Iplat.e, which were used to approxima- 
tely calculate the suction boundary layers by Schllchting [21], Head . 
[22], Eppler [23]. etc. 
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5.2 Integral Conditions 'for Momentum and Energy 

i 

, In this discussion, the rotationally-symmetric laminar and turbu- 

[lent boundary layers along the impermeable wall are calculated on the 

basis of the integral conditions for momentum and energy, with single 

parameter velocity pr.ofiles. The approximate theory for plane flows 
! } 

•:is discussed in detail; in [11], and therefore we will only give the 
'most important relationships here. 

I 
I 


1 I ^te.loc-ity.,profile of boundary layer 



x' 




; Using the boundary layer variables 
1 ■ ^ 
t 

^Displacement thickness 


momentum loss thickness 


■Energy loss thickness 
I 


Shape parameter 




(5.1) 


0^ 





L J 


Hi2 


cTi 

cT. 


; T. 


the- Pr-andtl boundary layer equation .is the following for stationary, 
incompressible, and rotationally-symmetric flow 




! B u , 3 u dp 


+ 


Sr 


dx' ■ 0y 

f 


r = 




9 u 
e 9y 


(5.2) 
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:and the continuity equation is 


B(ur)^ 3(vr) _ „ 

“57“ ■ ' 


(5.3) 


'with the boundary conditions 

y = 0 ; 

I ■ y = cT: 


u = V = 0 
u = uj-(x') 


jThej derived integral conditions for the momentum are 


2 ^ V 1 duj- , ^ 1 dr ^ 


(5.^) 


( and for the energy 


3 H. 3d-, -L ^ + <r, - - 2=n ' 0 

dx 3 Ujj- dx • 3 r dx D 


(5.50 


; where 


local friction coefficient 


dissipation coefficient 


- oI *T^du 
3“ 


C5.6) 


It is assumed here, that 6^/r |<<'.-l 5 that is^ at the nose tip and /17 
in the immediate vicinity 'of the rear tip^ the boundary layer conditions 
are violated. 1 


By introducing the new variable|s , we have] 


Thickness parameter 


shape parameter 


T' n 

Z = O Re cf _ . Re<p_ = 


|_| ^ ^ 
32 cTo 




C5.7) 


Then instead of 6 ^ and 6^ the integral conditions can be converted- 


aS' follows : 


) 



Momentum theorem: 


dZ 

dx' 


+ 


_ 1 _ 

u<r 



+ (n+1 ) — 
r 


dr 

dx' 


) Z - F, 


0 


C5. 


Energy Theorem 


d H32 p J_ duj- 

dx' 3 dx' 


H 


32 


- Ia = 
z 


(.5. 


with the abbreviations 

, = 2 + n +(n+l) 

F^ = 0+n) Re 2-2 

*"3 " ^"^12 

^4 ^ ^ 32 ^ ^2 ' 

/ 

If a single parameter velocity profile is used, the expressions 
C5.6) can be written as follows; 



"f 

"d 


Red 2 


B 



*1 




C5 


where the quantities a and 3 only depend on ^-^ 2 ^ which can be seen 
for a laminar boundary layer from the conversion 


oc 

B 


^2 

cT 

cT 


a(u/uj) 


Jw 


3 (u/ Ucf ) 




and for turbulent boundary layers. 
In this way, we obtain 


this is confirmed by measurem.ents 


F^ = 2 + n + (n + 1) H ^2 
F 2 = { 1 + n) CO 


*"3 ^ " *^12 

F4 = 2 B Re ^2 ■ ^32 


C5 


8 ) 


9) 


10 ) 


11 ) 


/18 


. 12 ) 
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Since for a single parameter \relocity profile, ^^.n also be 


represented as a function of ^^ 2 ^ 
considered as universal functions 


the abbreviations (5.-12) can be 
of E ^2 with the coefficients 


laminar: n = 1 

N = 1 


turbulent: n = 0.268 

N = N(H22) 


(5.13) 


! 


In order to solve the boundary laj^.er equations. (5.8) and (5.9), we 
therefore also need the relationslilps for the velocity profile 


12 


\ = f (H32) 
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5-3 Relationship's' for the VeTocit'y Profile 
5 . 3.1 Laminar Boundary Layers j 

No Suction . 

In the case of a laminar boundary layer without suction action, the 
velocity profiles are approximated! by the similar Hartree solutions of 


the type u^(x') ~ x 
m = m ( X ’ ) . 


m 


for the arbitrary external flow, and we set 


The Hartree profiles are defined in the range 


H 


32 


= 1.515 according to a = 0, m = -0.0904 (laminar sepa-r 
ration) 1 


H„., = 1 . 638 , according’- to m -> «> 
3^ ! 


and 


For the profile relationships*. required, we have the ¥alz [11] 
approximations ! 


06 = 1,441 (H32 “ 1,515)°'*^*^° 

(5 = 0,1573 + 1,691 (H32 - 1/515) 

= 4,030 - 4,183 ‘ 

L 


1,637 

0,3945 


(-5.14) 


which are shown in Figure 2 for the entire Hartree range. 
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; The more accurate relationships can be deterihined from the single 
j parameter formula of Geropp [l6] for the laminar velocity range^ as 
■ shown in Figure 2. | 

; With suction | 

! ^ 1 

When the laminar boundary layer Is influenced by suction^ the simi- 

1 

• lar solutions for the plane plate [with homogeneous .suction are used for 
! the velocity profiles, which differ substantially from the Hartree so- 

t lutions. In. order to have a good | approximation in the entire' range, 
! Epper [23] approximated the Hartree profiles with a pressure increase, 
and the sunction ,pro.file.s..,of-^-;Iglisch ,[.19-]: In .the .case of a pressure 

* drop, '..by. using the limiting case of asymptotic suction profiles = 

' 5/3s as shown in Figure 2. j 


oL = 2,512589 -1,686095 ~ *^^2 

1,51509 ^ H22^ 1,57258 

c<, = 1,372391- 4,226253 H + 2,221687 


^5^1,57258 < H, 


1,66667 


6 = 7,853976 - 10,260551 H32 + 3,418898 H32 

^ for I 1,51509 ^ 1,66667 

4,02922 - (583,60182 - 724,55916 H22 
+ 227,18220 H22) 
fT5^ 1,51509 < H32^ 1,57258 

H^2= 79,870845 - 89,582142 + 25,715786 H22 

[f3Fl 1,57258 < H^2^ 1,66667 


(5.15) 


The separation profile (ct = 0) is specified by -l,51509jfr^ 

j ' ' . 

Using the approximations (5.15), very reliable results are achieved, 

' I 

even for a pressure increase. If one calculates the delay stagnation' 
point flow 
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Uf ^ 1 X 

u 1 

_ qp ^ - 

for which there is an exact solution of Howarth [24] available, the 
parameter separation 'poirTt to an accuracy of Q.2%, as shown in Figure 

3. 

In the case of continuous suction with arbitrary suction inten- 
sity, no completely -laminar velocity profile can be obtained as 
5/3, behind a suction slot there are velocity profiles with R ^2 ^ 5/3 
as well, which violates the Prandtl boundary layer assumption v/u << 1. 
With the exception of the R^^ approximation, these • relationships (5.15) 
can also be used for R^^ ^ 5/3-' The approximation /21 

H^2 = 3,738 - ■y 13,43 . (5-16) 

for Hj 2 > 1,6667 

contains the rectangular profile with R ^2 ~ 2, = 1 as the limiting 

ease . 


5.3.2 Turbulent' Boundary Layer 

According to Fernholz [l8], in' the ease -of a turbulent boundary 
layer, the following semi— empirical laws can be used with success for 
the velocity profile over the Reynolds number range of interest: 

Wall shear stress law according' to Walz [11], 

06= 0,0566 H 22 “ 0/0842 (5.17) 


Dissipation law of Rotta [25] and Truckenbrodt [26] 

_ ‘ 

13 = 0,0056 
N= 0,168 

i 

Shape parameter relationship of Fernholz [18] 

H^2= 1 + 1,48 (2-H32) + 104 (2~H32)‘^'^ 


(5.18) 


(5.19) 
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Separation occurs for E ^2 


= 1.50 - 1.57. 


;5.^ Stepwise Solution of Integral! Conditions 

' I 

Using the relationships givenjin the previous chapter for the 
velocity profile^ the boundary layer calculation is reduced , to the 
Isolution of the equation ’ system (5i8) and (5-9), which applies .both 
for the laminar and turbulent boundary layer, and the two unknowns 
^are Z(x’) and H 22 CX’). j 

I Assuming that- in a small interval we have 

- AxJ ^ = xj - X._.| 


;and that the body contour and the velocity profile can be approxi- /22 
-mated by a linear law _ 1 . - - - -- — 


r. “ r. , 

— r , . "t- “K 7 (x "" X, . ) 

i-i Ax. . , i-l 

i-l 


U ^( X ') = + 


Vi “ Vi-l 
Ax: . , 


( 5 . 20 ) 


(x' - x.;_^ ) 


iand if we introduce average values jfor the variables ZCx-') and H„.,Cx') 

1 0-^ 


land for the universal function P (x'), 

I .y V / J 


■Zi,M 

; ^325, i-l " 2 ^*^32i ^32i-l^ 

i, i-l " ^^32i, i-1^ ' ^ " 1, 2, 3, 4 


C5.21) 


jthe momentum and energy theorems can be solved in closed form. The 
'step formulas are the following in jdimens ionless not at ion in the ro- 
jtationally symmetric case I 


r. ./I ln+1 Z. , 

I-l I-l 


_T - Az - 


*^32i ’^321-1 


Ax: . - 
1 , 1-1 


Ax: . VI 
1/ I-l 


r. ,/I 

1+ ' y. 1n+l 


■C5.22) 


■with the coefficients 
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= 


A.. = 


H 


= 



u^. ,/u 

K 1-A^ 

- — ^ — L.y.<j/^CQ 

’ 1 - ■ 


C5.23) 


1-A 


i/u 

d 1-1'' oo 


i B 
1 H 




rTK 


1 ' CO 


1 - 


U r. ,/u 

° i-r GO 
Ux./u 


oo 




and the step along the' contour is 


^\x. , . - 1 1 X. X, - 2 r. r. . 2 

- ll(^- _l^) + (t^- -fi) 


C5.24) 


1 ' r._,/l ■ ■ 

In the plane case, in (5.22) we have. . ! y, ■ — —I,, and In (5.20) and (5.24), 


is replaced by r^^. t 

The stepwise determination of 
H ^2 at the end of the. interval (i) 


71 


the boundary layer parameters Z and 
is done by an iteration from the 
values' at the beginning of the interval (i-1). Betails are given in 
Chapter 8.2.2. j 

The boundary layer calculation for a laminar boundary layer is 
started at the following points : | 

- nose tip of body of revolution 

- suction slit 

and for a turbulent boundary layer, 

i 

- the laminar-turbulent transition point 

The end of the boundary layer calculation occurs at the following 
points for the laminar boundary layer ' 

- one suction slit ; 

— — the- l-ami-na-r.-turb.ul.ent^=tr.ansi..tj.-on. . point. 
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, - or the laminar separation point. 

■and for turbulent boundary layer, 

< 

- the end of the range of validity of the integral conditions in 
the area near the tail tip of the body of revolution 
I - or the turbulent separation |point . 

I 


i 


5-5 Initial Calculation I 

; At the nose tip (x = 0), the step formulas G5-.22) fail because the 

' . t 

'thickness parameter is Zq/1 =0. ] 

Therefore, in“"the' f irsi'etep 'Ax"' ‘ ^yi the step - formulas C5-22) at 
I the nose tip are replaced by relationships which are derived from 


conical flow: 


J. 


^ 32 , ^32 ^32 

I o 


]_ 

3 


l+rrig- ”1 


C5.25) 


j and the cone angle and the Hartree iparaitieter of the plane replacement 


flow 


I 


m. 


1^ 




C5.26) 


are specified by the cone angle 




— arctan — — tt 
ir • x^/1 


C5.27) 


In the plane case, the relationships for the first step are: 




^32 ^32 ^32 

1 o 


F x' 
r-2 


( 5 . 28 ) 


where 


l+mF| ~r 


% 




2 y/' 

- — arc Tan tt 

ir / 


(5.29:) 
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mTFii the assumption that Hartree velocity profiles occur at the nose /25 
tip, the relationship and 3') can be approximated by 

= 1,5720+ 0,1561 6^ - 0,2470 1^ + 0,2244 - 0,0804 

^ b t b b 

for^ 0 < 1 (5-30) 

> 

H^2 = l,5720+~0,2258 ^^ + 0,7663 + 5,3850 

~for~ -0,1988 < < 0 


The stagnation point flow 

rotationally symmetric (6=1; fi^ = 0,5) : H 22 =1,6113 

plane (3 = 1) 13^2 =1,6250 

and the plane plate flow 

B = 0 : H32^ = 1,5720 
are contained asl special cases. 


5 . 6 Determination of Friction Drag and Pressure Drag 

Because there is no boundary layer separation along the tail part, 
the external drag that is the sum of the friction drag and pressure 
drag of a body of revolution can be determined from the momentum loss 
in the wake. According to Young [27] and Pretsch [28], the drag coef- 
ficient of a body of revolution is the following with respect to sur- 
face area 



Wn 2 Qqo 

qco o ~ o 


(5.31) 


with the following definition of momentum loss area 

© = 2'tt f — (1“- — ) ‘■'^'' (5. 32 ) 

J uj uj- 

r 

By using the integral conditions for the momentum C5-^) on the wake 

(t^ = 0), the momentum loss area 0^ at a large distance behind the body/26 
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(p = Poq) can be reduced approximately to the boundary layer variables 
,to the tail tip 


5+H 




oo 




oo 


12h 


(5.33) 


In the immediate vicinity of t-he tail tip, however, the boundary 
layer calculation fails because of |the integral conditions (5.8) and 
(5-9). This is because the condition 6/r << 1, is not satisfied. It 
*is found that the boundary layer calculation must be terminated (x = Xg) 
iin the case where ^ 1/15. 

- In this investigatl'on“ 

I - laminar boundary layer suction 

i 

I - high Reynolds numbers 

^ - thick bodies 

'however, we have 1 — -x,.,/! << ' 1 so that the friction drag component of 
. ii ' 

;this region - 


A 


AW, 


c = 

Wb 


R _ 


J 

4tt 


r 


R O/! 

L_. _ _%/' 


(“) 


oo 


r 

r 


Cf d(]i) 


due to 


'can be ignored, compared with the .t^otal drag coefficient 
*the small local wall friction and the small body radii in t*he vicinity 
of the tail tip. In addition, in the case of bodies with sharp tail 
tips, the shape parameter in the tall tip region is almost constant 


= H 


= H = const, for 
12h 12 e 


0 . 


fTheh! according to (5.33)j we have i 






5+H 


12e 


C5.3^) 


,and consequently. 


0 


oo 


5+H 


12e 


= ©u (rr^) 


CO 


C5.35) 
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so that the external drag coefficient is 




w. 


= 2 i^) 


5+Hi2, 


O ' u 


(5.36) 


oo 


According to Scholz, [29], we have the following relationship for 
the moment 11111 loss area 


=2. ^ik (, 4 ^) 


~T 


I 


2 


(5-37) 


so that the drag coefficient can be expressed by the boundary layer 

MI2e 


variable at the position x^/l 


4 'TT 


cP? 


w. 


o/r 


where 


^■2. 


, (C /I 


(^^ReY 


(5-38) 


L U 


n+1 


oo 


^12e ~ ^^^32e^ according to (5.19). 

In the plane case^ .the boundary layer is calculated to the trail- 
ing edge. The drag coefficient expressed by the trailing edge varia- 
bles is then given by 


where 




2 


“T 



oo 


5+Hi2h 

2 



Zh/I 

Re,)" 

u I 

00 


1 

n+1 


^12h " ^^'^32h^ according to C5..19) 


C5.39) 
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' The wake behind a flow body is turbulent at high Reynolds number^j 
iRven. 1 for laminar boundary layer suction to the end of the body, one must 
lestablish the latest boundary, layer transition for the drag calculation. 

, Body of revolution: ■ iiji- <• 


plane contour: 


< 1 


Comparison calculations for the drag determination are given in 


'Chapter 6.5. 


6. Determination of Instability Point and Transition Point Over Bodies 
of Revolution. i 


;6.1 Preliminary Remarks 1 

' i 

The laminar-turbulent transition point over a body in the flow 
j could not be exactly calculated up po now, because the way in which 
I turbulence is produced is not clarified. 

[ In experiments, Schubauer and tSkramstad [30] showed that [boundary 

I layer transition must be attributed' to an instability in the laminar 

i boundary -layer . | 

' ' ! . . 

; If the disturbances in the laminar boundary layer -caused by rnci- 

jdent flow turbulence sound, vibration, or surface roughness are small , 

t 

Sthen the transition process can be 'divided into three regions: 

^ 1 

, I. stable laminar boundary layer 

, All perturbations are small in amplitude, and affect the lami- 

' ■ • ! 

; nar boundary layer and decay in time, 

, II. Unstable laminar boundary ilayer. 

' At least a few partial oscillations of the perturbation motion 

t j 

are built up. In the terminal phase, secondary and high fre-, 

quenc.y perturbations are superimposed, which forms -a three- 

dimensional vortex system.j 

-- III-. Laminar-turbulent transition range 

: Due to an instability in the secondary flow, turbulent ■ spots’ 

• are produced in the laminar flow, which propagate during the 

, • ! . /.2 
downstream motion until -the completely turbulent state is 

i reached. I 
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.! The determination of ’the instahi-lity point after which there is a 

) buildup of the disturbance motion p diseussed in stability ■ theory . 
i Mathematically .the problem consists of solving the Orr-Sommerfeld per- 
/turbation differential equation, which, was successfully done, first by 
[Tollmein [31] and later on ^ Linj [32] for the Blasius profile .of the 
, ^Tan^ boundary layer. VJhen there is a changing pressure gradient, the 
' instability of the laminar boundary layer was calculated by Schlichting 
i and Ulrich [33] for the Pohlhausenj P6-profile . This was done by Pretsch 


[34] and Tetervin ['35] for Hartree 


profiles 


The stability of the laminar boundary layer with suction was in- 
vestigated by Bussma-nn--and -Muenz [36'] for -the as-ymp-totic suet ion... pro- 
file, It was calculated by Ulrichl [37, 38] for the -similar suction 
profiles given by Iglisch and by Schlichting and Bussmann. As long as 
the transition process is not completely clarified, only empirical re- 
lationships can be given for the Ipgth of the buildup length. The 
transition region, however, is so port according to experience that 
in the first, approximation it can be ignored, and can be replaced by 
1 the transition point which refers to the beginning of. transition. 

I Since in the ease of small per'burbations , the transition point 

! lies between the. instability point! and the laminar separation point, 

! these two limiting positions are open used as the transition condi- 
; tions. In the case of high Reynolds numbers, the distance in the 
j instability point and separation point is very large. 

' More recent t-ransition criteria are based on stability theory, 

i Michael [39] gave an empirical transition criterion for plane jboundaryj 
! layers which later on. was analyzed] theoretically by Smith and Gamberoni 
j [40] using the build-up, diagrams of Pretsch. [4l]. The criterion of 
> Granville [42] also considers the influence of the pressure gradient 
' on the length of the buildup path. In the case of plane flow, it is 
: possible to predict with some success the 'transition point using these 
f criteria, but information for bodies of revolution is not reliable. 


6 . 2 - 


-Stability Criteria ' 

' Pretsch [43] showed that the results of the stability theory can 

I be transferred from the plane casej to the rotatlonally-symmetric case. 

as long as 6/r << 1. | 

' Figure 4 shows the calculated] instabilities 
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- for the Blasius profile of the Ip lane ' boundary layer of Tollmien 

[31] and Lin [32], ' 

- for the profiles of a boundary layer with a pressure gradient of 

* Schlichting, Ulrich [33] for Pohlhausen-P6 profiles 

* Pretsch [3^] and Te tervin [35] for Hartree-profiles 


in the form 





The various stability calculations ' agree well. The pressure drop 
(H 22 > 1 . 5726 ) has a stabilizing effect on the laminar boundary layer, 
whereas a pressure increase ^ X. 572.6) has a deep stabilizing ef- 

fect . 

With suction there is an even greater stabilization of the laminar 
boundary layer. 

The results of the stability calculations of 

- Bussmann, Muenz ■[3&‘] of the asymptotic suction profile 

- Ulrich [ 37 ] j for the Iglisch suction profiles, 

- Ulrich [ 38 ], for Schiichting-Bussmann suction profiles, 

are also shown in Figure 4. 

One sees that the influence of a pressure gradient and suction on 
the stability of the laminar boundary layer can be approximately re- /31 
placed by a single criterion. The following are used for approximating 
the stability criterion: 

- in the case of a pressure increase, the results of Pretsch ’[3^] 
for Hartree-profiles, 

- for a pressure drop, the results of Ulrich [37] for Iglisch 
suction profiles with a limiting case of asymptotic suction profiles. 

Approximately, one obtains the following for the position of the 
instability point as a function of ,the shaped parameter 


7 ^ /f 0 

'log Rerf-2^ = 4,556 - 76,87 (1,670 -H32) ' (6.1) 

in the range T,'5l5'o '<Tl~2‘'^'“r,6M77” 

Using the approximation (6.1)',, in the case of a pressure drop, 
the influence of a pressure gradient is very accurately represented 
up to the r’otatlo.nally-symmetric stagnation point CH ^2 = 1.6113). 
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6. 3 Transition Criteria for Plane jPlow 

, 6 . 3.1 Discussion of the Known Criteria 
! I 

' By evaluating the transition measurements which have become known, 

.empirical transition criteria havejbeen found which are similar to the 

i parameters of stability theory. It is assumed that turbulenc/e is pro-,/ 

! . I 

duced by the fact that the perturbation amplitude reaches a certain 
‘level. The buildup of the perturbations '.depends considerably on the 
'pressure gradient of the flow. 

! Granville [42] assumes that for plane flow, the length of the 

; buildup path characterized by 'the Reynolds number difference 


•|ARe^2 ~ 


jonly depends on the average pressure gradient betvjeen the instability 
j point and the pressure point, withjthe assumptions made. It is given /32 
by the following expressed in terms of the average value of the Pohl- 


I hausen shape-parameters. 


.A. 




I 

I where 




K 


Xy-X^ 


'Vc(x') dx' ; 


( 6 . 2 ) 


3^(u/uj) 


w 


Prom the corresponding evaluated transition measurements, the em- 
Ipirical law • 1 . ' 




C 6 . 3 ) 


I 

I 

jhas been derived. 

, Smith and Gamberoni [40] assume., that the self-excitation of the • 

*Tollmien-Schlichting-¥eilen waves re^^esen^ the dominant process be- 
' tween the instability point and the transition point, even though the 
assumption of small perturbations is no longer satisfied in the vici- 
'nity of the transition point. They introduce an apparent amplification 
ratio, at the beginning of transition, fto which the anplifcudes of the 
Tollmien-Schlichting waves are to ije subjected with the prevailing 
critical frequency. Using the built-up diagrams established by Pretsch 
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■ [4l] for Hartree profiles in [40] , jnumerous transition measurements 
■ have been evaluated. It has been found that the apparent, ampllflca- 
Ition ratio | 

; n = =8103; (c. 


a = e‘3 


8103; 


(6.4) 


can be used as the transition parameter, where 


the damping buildup is 


< 0 damping 


’O” ’buildup' 


^ This finding was also confirmed by Ingen [44]'. 

With this convention, it was possible to subsequently correct /33 
tthe[pur^I^ empirical transition criterion of Michel [393 in [40] 


The criterion was approximated as follows in [40J 


0,46 

■■ i JT 7 

for 3 • 10^ < Re <2-10 


(6.5) 


The bodies of revolution, the transition criteria C'6.3) and (6.5)', 
are no.t suitable,- because the buildup ratio is expressed by boundary 
layer. variables which are different in the plane and rotationally-sym- 
imetric cases. The. suggestion of Granville [42] to use the criterion 
'(6.3) for plane boundary layers injan unchanged form for rotationally^ 
’symmetric boundary layers must of necessity lead to inaccurate results. 
(However, Smith and Gamberoni [40] point out that in the case of bodies 
'of revolution, the criterion (6.5) must be applied to the plane replace- 
ment boundary layer. However, this does not consider the fact that the. 
boundary layer transformation can only be performed along the buildup 
'path. j 
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Because of the deficiencies of the transition conditions mentioned 

t 

above, we suggested a semi-empirical criterion which can be transferred 
from the plane case to the rotationally-symmetric case. Since the tran- 
sition points will be less reliable, the more single processes of the 
transition process are considened,! and the instability point can be de- 
termined very accurately, we will ^consider the buildup path separately. 
Therefore, in the (following /, instejad of considering the transition cri- 

t 

terion, we will give for the length of the buildup path. 

Using the' transition condition (6.i0 we can derive a criterion for 

the length of the-buildup pa-th-in -hhe' case' o-f similar velocity profiles 
/\ I 

C6 or H 22 = const.) between the instability and the transition point /3^ 
using the buildup curves of -Pretsch [4l]. In this, the influence of 
the pressure gradient and the suction [isU included. 

According to [4l]., Figure 5 shows the following buildup curves 


•r= 0: stability condition according to [34]; see 
' 8,<jb ■ ■ “ 

■ ^ ■[ = 9: transition condition according to [40] and 

- [4#]lj-- ■ 

I 

For a constant-shaped parameter, the distance between both curves 
expressed by the boundary layer variables given in Chapter 5,. where 




( 6 . 6 ) 


where , = f{B) = const, according to ( 5 - 30 ) 

. ' 


t 

is a measure for the length of t he j buildup path. 

The evaluation of the transition measurement for wing profiles 
with low turbulence initial flow [30, 46-54] done by Moeller L45] shows 


that the relationship (6.6) ^s,/a mere/ reliable criterion for the length, 
of the buildup path for constant shape parameter than the one discussed 
in^6.3'vl, if we introduce the follLwing average value for the shape 
parameter for in stabilit y an d t rans ition 




^2j 


[([432) 




H 


32 xyi - x^TI" 


H32 (f)d(^) 




C6,7) 
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This is shown in Figure 6 . The criterion has been confirmed by measure- 

I 

ments if one ignores certain scatter, which is due to the differences 
in test conditions of the individual authors. 

With special consideration of the transition measurements [30], 
criterion ( 6 . 7 ) is approximated by the following expressions: /35 

log (Rej- 2 ^ - = 1,6435 - 24,20 {] , 5]50 

for”" 1,5150 < ^22 ~ 

( ( 6 . 8 ) 

— 2 715 

, log (Re,T 2 ^j " " 3,312 -967,5 (1,6250 -H 32 ) 

for~ 1,5600 < H.o ^ 1/6250 - ' 

This is shown in Figure 6 , so that the transition point of a plane 
laminar boundary layer can be determined directly. 

The criterion of Granville [42] shown for comparison in Figure 
6 only gives approximately the same values in the central shape para- 
meter range. 


6 . 4 Transfer of tlie' Transition Cri'terioh to Bodies of Revolution 

Smith and Gamberonl [40] determined the the transition cohdlton 
(6.4) also applies for rotatlonally-symmetric boundary-layers. This 
means that criterion (6.7) or its approximation (6.8) can be trans- 
ferred to the plane replacement boundary layer of a body of revolu- 
tion using the Mangier transformation [12]. 

Using the assumption 


( 6 . 9 ) 

( 6 . 10 ) 


the coordinates of the plane replacement flow are the following for 
a body of revolution with the contour r(x'). 


y 



(x') dx' 




y 


( 6 . 11 ) 
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where the constant L is a characteristic length. However, if we intro- 
duce L = L(x') in such a way that everywhere • 736 


x' = x' 


( 6 . 12 ) 


then we obtain 


dx' 


(6.13) 


The boundary layer transformation starts at the instability points, so 
that 


and because of (6.11) 






Consequently, the Mangier constant is calculated from the following 
using a dimensionless notation 


f ifi d(f) 


(6.14) 


in the 


x' ^ x(, 


range < -p < 


Using the assumption (6.10^ and the relationships 

'J' - cT 
"^3 ~ i/r 

according to (6.11), the transformed thickness parameters are 


fed-j =■ Rejj ; Refj = -f h 


(6.15) 


and the transformed shape parameters are 


^32 - *^32 


^32 ^32 xVl - x"! 


/ ^ \ I / -A ' 


H32 


(6.16) 


32 



2 I 

The expressions (L/I) according to (6.l4) and ^^2 ^-^cording to /37 
( 6 . 16 ) are formally the same, and are calculated in steps according to 
the boundary layer parameters (5.7). This means we have the following 
step formulas in dimensio nless .1 form for the Mangier constant 



(6.17) 


and the average shape parameter is 


^ 32 . ^ 32 . , 

I 1-1 


Ax'. . . 

',■-1 

I 

Ax'. . , 
1,1-1 


H + H » 
32j-l 32; 

2 


The step along tiie contour AxJ . -./I is specified by (5.24). The 

-L ^ J_ X 

instability point is the initial point. 

The transition measurements for two bodies of revolution with a 
low turbulence Incident flow [54 and 55] have been evaluated using the 
above method, and this is shown in Figure 6. Because of the fact that 
the measur.ed points of the bodies of revolution do not have a greater 
scatter than those of the profiles, we therefore have found an experi- 
mental verification of the transfer rule for the criterion (6.7) pf two, 
bodies of revolution. 


6 . 5 Comparison Calculations 

The comparison calculations were performed using the computer 
program discussed in Chapter 8. 

The criterion for the instability point and the transition point 

i 

can most easily be tested using a plate in longitudinal flow, as shown 
in Figure 7 . The stability criterion (6,1) contains the exact solutions 
for the Blasius profile of Tollmieri [ 3 I] and Lin [32] in a satisfactory' /38 
way. The measurements of Schubauer _andl Skramstad [30] are well-repre- 
sented by the transition criteria (6.8). The calculation drag coeffi- 
cients for natural transition and completely turbulent states agree 
well with the laws given in [10]. ' 

The transition measurements and the drag measurements of Boltz, 
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; Kenyon j and Allen [5^s 55] are good test examples -for the rotatlonally- 

symmetric case. As Figure 8 showsl, there is good agreement between the 

! calculated results and the measurements. The theory predicts somewhat 

'longer buildup laminar paths at' hi^h Reynolds numbers, and accordingly 

I at lower drag coefficients. A calculation using the calculated transi- 
' -I 

' tion point positions showed that there was no difference to be -detefc ^ 
, ted between the calculated and measured drag coefficients. 


7'. Influencing t lie Laminar .Boundary .Area by .me.ans of Suction Slits 


7.1 Preliminary Remarks | 

Suction of the boundary layer! particles near the wall causes the 
following : j 

- reduction of the boundary layer thickness 

- the velocity profile becomes fuller, and 

- the potential flow is influenced by a sink effect. 

J 

By changing the boundary layer variables, a stabilization of the 
laminar boundary layer is brought about, as can be seen from the sta^ 
bllity diagram in Figure h . t 

Even though a continuous distribution of suction over the sur- 
face jwould’b^f optimum from a suction power point of view, only a local 
suction through slits arranged perpendicularly to the flow direction is 
possible in technical applications-. Also gaps or perforated, strips 
could be- used, as shown in Figure 9. Since the influence on the lami- /39 
nar boundary layer is very similar| in all three eases, we will ’only dis~ 

ft 

cuss suction slits in the following. 

In the boundary layer calculation, the suction slits represent 
discontinuity points, because the assumption' of the Prandtl boundary 
layer equation at the wallv^/u^^ is violated because of the con- 

siderably increased suction veloci'ties. Influencing of the boundary 
layer at the suction points must therefore be treated separately. Two 
approaches are possible: | 

a) the suction points are considered as strips with a finite 
width, with an intensive area of suction. 

For this case, exact solutions have been given by Rheinboldt 
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'[56], Smith and Clutter [57] and Krause [58]ji as well as the approxi- 

;mate solutions of Bethel [59] and Wuest [60] for specific applications. 

b) the suction points are considered as slits w.ith a negligible 

Vldth. In this case, ¥alz [6l] and Wuest [62] suggest | to>| appr.oxlmate- 

,ly determine the boundary layer parameters behind the suction point 

>from the started residual profile s^amputation principle). 

The solutions- of a continuous j suction over a path-are complicated. 

‘They are also restricted to a small suction velocity range which is not 

lof practical interest.- The amputation principle can be used for an ar- 

'bitrary suction velocity, body shape, and pressure distribution. 

' i 

In the following. we^ will -gi-ve jthe- -gener-al-appro-ximate relationships 

^for influencing the boundary layer ‘.by means of suction slits according’ 

■to the amputation principle. These are valid both for the laminar and 


,the turbulent boundary layers. After this, we will apply this princi- 
!ple to the laminar boundary layer and will investigate the effects of 
jSUCtion slits on the laminar bounda-ry layer. 


I We will show that the sink effect of the suction slits is in gene- 
■ral negligible because of the low sucked-away amount, when the' boundary 

t I 

jlayer is influenced. _ ^1. . . ■ /40 


I t 

■ 7.2 Approximate Solution According- to Amputation Principle 



The amputation principle of Wadz [6l] is based on the following 
■idea: we will assume that a velocity profile I exists ahead of the suc- 
tion slit, which is divided up into[ the following: 
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slit; 


- a region' near the wall of height which is sucked. into the 


- an external region (residual profile) which goes along the stag- 


nation line, vrhich goes on through 


the slit. It is deformed without 


losses in such a way that the incident velocity u(y„) decreases at the 

{ y 

rear slit edge to -u = 0 .(wall condition) . It forms the initial profile 
II for the boundary layer development behind the slit. 


’7*2. 1 Changed Boundary Layer Variables Behind the Suction. Slit 

The changed boundary layer variables II behind the slit _can be- de- 
termined from the .residual. pr,of.lle.^-wit.h.--the-.assump-tion that the mass, 

! 


.momentum, and energy of the initiated residual profile do not change, 
I \ 

'during the motion over the slit, | 


If we select a single parameter velocity profile 


/11 




(With the profiles according to 




I _ 


‘i'l 

<h. 


I _ 


d; 




cT. 


3i _ 


cTr 




, and the initiation point is at the .relative height ■ ~ ^ then 

the boundary layer variables of the residual profile are 




Sr 


s 


^3i 

cT. 




a*. 

1 




(7.1) 
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In this way we obtain the changed jyarlabies for the 'boundary layer 
calculation behind the suction slit 


■II 


^32„- 


^31 / 


C7.2) 


In order to again order the velocity profile II into the family of 
the single-parameter velocity profiles, instead of the shape parameter 
H^ 2 ii will introdjic.e_ the .average d_,v^al.ue /42 


H 


32 


E 


9 ^*^ 32 - 


JL 


32, 


(7.3) 


where 


! 




h) 




C7.4) 


The functional relationship in C7-4) is unique for single-parameter 
velocity profiles. The. shape parameters H 32 J-J ^32n only' 

slightly from one another. 1 

In the case of single ‘parameter velocity profiles ahead of the 
suction slit, the boundary layer variables behind the slit only depend 
onthe shape parameter thet relative suction height tIq. 


*7,2.2 Suction 'Variables | 

' By using the following definl|:lon of the relative suction thick- 


ness 


cTq _ f /U \ J = n - 
0 


C7.5) 


we can calculate the amount of alrj sucked away by a slit corresponding 
to the suction height referred! to the slit length b (body of revo- 

y I 

lution: b = 2 irr)'. ' 
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Q 

b 




or the local suction coefficient amount 


Q 

. u B T 


CD 


^<r *^0 
u 

OD 


(7.6) 


In addition, we define the nondimensional suction amount coefficient 

Re. 


'I _ AZi 


Ref 2 <f2i A 


where 


(7.7) 
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Re, 


Re 


Re 




I 


which only depends on ^ single parameter velocity pro- 

file ahead of the slit. 

According to Gregory [63], the disturbance of the boundary layer 
due to the slits is minimal if the slit s is equal to the suction height 


n 


Q' 


Accordingly, the optimum slit width is 


T = ^ “T 


C7.8) 


The,' nondimensional slit width 


s Re I 




T Rerf-2^ 




C7.9) 


only depends on single parameter velocity profiles. 

Using the continuity law, we obtain the average suction velocity 
in the slit as follows: 



(7. IQ) 
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7 . 3 The Amputation Princlpl'e to the Laminar Boundary Layer 
7-3.1 Suction Relationships 


The relationships in Chapter 7-2 apply for any boundary layer 
state. Here, we will derive relationships for Influencing the lami- 
nar boundary layer by suction slits. ¥e will assume the laminar velo- /44 
city profile ahead of the suction slit can be approximated by the single 
•parameter model of Geropp [1&] ' ~ ~ i 


= l-(l-r^)^ (l+a^r^ + -f 


C7-11) 


with the coefficients 


= c - e 
_ A 


2 

_ A 


- c s + 


(c+1) 


^3 2 2 6 

and the following analytic expressions for 


= - 0,89855 + 111,60901 A + 13,6795 


^ r B (u/ u^r) 

L Jn." 

= 5,22550 + 1/1,30839 A + 10,85171 


c (c+1) c (c+1) (c+2) 
C - £ + 


(7.12) 


(7.13) 


The Pohlhausen-shape parameter which occurs in these relationships is 
related to the shape parameter 


A = 


9 ) 



I'oy means of the following approximations 


A= 38,745 (3 - 7,1178 6^ + 6,3726 


(7.14) 


and 


O Q 

8 = -2788,62 + 5439,97 H 22 " 3539,62 + 768,15 

for 1,5237 ^ ^ 1,5729 

= -35461,01 + 67123,15 H 22 " 42356,63 + 8910,45 

for 1,5729 < S 1,6239 


I (7.15) 
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The Geropp model is valid in the following range: 


-0.1988 5 B < 1.0 

Separation profile 

1.5237 1 1*6239 


> i plane i 

, [stagnati o nf 
! point profile-./ 


The boundary layer thickness used as the reference yariable for 
the suction height is defined in many ways*: 

- in the boundary layer theory of Prandtl = y(u/u^^ = 0,99) 

- in the approximate model of Geropp = y(u/u^ = 1,0), 

One must distinguish between the following relative suction height 





of practical interest, and the theoretical quantity 


Iq cTe 


which occurs in the suction relationships. 

The relationship q 

is shown in Figure 10 and can be approximated by the following polynomial 


~ = 5378,67 12707,78 11189,97 - 4346,14 + 627,10 

, cT'^ ' 32 32 32 ( 7 .I 6 ) 

for 1,5237 < < 1,6239 . 

The boundary layer variables ahead of the suction slit are written 
using the -velocity model of Geropp 


*The quantity 62 is introduced as a reference quantity for the suction 
relation ships independent of the selected velocity model, but the eval-. 
'nation . of the following integral expressions would become much more 
difficult . 
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The integral expressions 

= f’ 

0 

I 


=1 (1-r^)^ (1+a^r^ + Q2r\^ + dr| 


^2 1 


I 


3i 


(1- f|)^^(l+a^r( dr| 

(1- r|)^^(l+a^f| + 




C7 


are evaluated in the appendix. 

Accordingly 3 we convert the boundary layer variables of the re- 
maining profiles 


%= = 1. 






E 


cf. 


3 = [ i^0-^)dn = I. -I. 

J uj- I Ijl 2]i 


cT- 




3 e _ 


Sr r< . 


d ^ 21 -31 +I 

J Uj- L Urf- J I In 2n 


3n , 


Q _ 


cT, 


'JcT 


d n = 


\ - <Ilx - ^1. ) 


the boundary layer parameters of the remaining profile 




to 

!| 


1— 1 

1 




C7 


C7 


.17) 


.18) 


.19) 


. 20 ) 
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H 


/ - *^311 /^i ^ 2 - ~ 

32]i 6'n 

2n T 


\-h^ 


»h, = 


12 


ii /<^i 


I 


■f2i/ 'i'- 


ll 4 

I ‘E ~ 


u ) 


and the nondlraensional suction variables 

6 


0 


_ aa + ^ln-^li 


'Iq 

75757 


iii - b, 

G 

III - l2i 


The relationship! = f(H ^ ) in (7.20) is approximated 

3/h U 5 


H'' = 2,660594 - 0,922220 H,„ + 0,276058 H,, 

32u o i 1^: 

- 0,036478 H,„ + 0,001836 H,„ 

IZu IZji 


,for; 1,5237 ^ H" ^ 1,5726 

dZji 


H" - 1,741773 - ^0,038887 H.. - 0,072234* 

32j[ ’ \Zji 

i for'; 1,5726 < H"- ^ 1,6667 analog (5.15) 


H" = 1,4418 + 0,07444 (3,738 - H,. T 

32 j£ I XI 

fori H"„ > ] ,6667 analog (5.16) 

32jj; 


The integral expressions 




b, ° J ('+“,T^ '■2^ + “s'l''*" 


3,2 


-E n; 


= j'(l-fl)^^ (l+a^t| + a2i]^ + a2r^Y d 

1ft 


3,3 


= ^32/ fa) 


are evaluated in the appendix. 


(7 = 20 ) 
(continued) 


(7.21) 

by 


C7.22) 


( 7 . 23 ) 
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7-3-2 Results Tor 'Arbitrary P'ressWe Gradients 

Figure 11 -shows the velocity jproflles behind a suction slit which 
was produced for different relative suction heights from the Blasius 
profile ahead of the suction slit.! Much more full' velocity profiles 
occur already for small suction heights. 

For arbitrary pressure gradients, we have plotted the following; 

- Figure -12 5 the shape parameter behind the suction slit 

- Figure 13, the ratio of the, thickness parameter before the 

suction slit, | 

I 

- Figure l4, the 'nondlmensional suction' amount coefficient c 

- Figure 15 j' the optimum nondimens ional suction’ width, s*. 

As a function of the shape parameter ahead of the suction slit H 
and the relative suction height ' , the shape parameter H 
in the limits 


32i 


Q 


32i 

can vary 


t 


t 


t 


1,5237 ^ H 


32 1 


1,6239 


for which Geropp velocity- model is! valid. 


suction removal 


3q 


= 1 . 


For complete boundary layer 

= 2 ) 


one obtains a rectangular profile 

— ' j - -- 


II 


with ^ 0 as a limiting-lvalue behind the suction slit. In 

addition. Figures 12 to 15 show the asymptotic suction profile as a 

f 

limiting case for continuous suction. 

In the case of the plate boundary layer = 1*5726)', there is 

a relationship between the suction' amount coefficient cf and the sue- 

I y 

tion height • which is shown in* FigiJire l6. Since in this case there 

> Q I 

is only one free parameter", the shape parameter behind the suction slit, 
the thickness parameter ratio and the slit width s* can be repre- 

sented as functions of the suction^ amount coefficient c„* as shown in 

1 - Q 

Figure 17 . 


In order to evaluate the slit 


suction, it is important -how fast 
The variation of the 
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the suction effect drops off behind the slit. 

1 , 

boundary layer parameters 2/1 behind a suction slit shown for 

a flat-=plate as a function of -the suction amount coefficient c shows 

y 


the following; 

- The shape parameter increase decays behind the suction slit, as 
the thickness parameter reduction still has a large effect at 'a large 
distance behind the slit. 
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- as the suction amount increases, the shape parameter increase 
^along the trailing slit .edge increases, but it decays faster because 
of the related thickness reduction.- 

t I 

This is shown in Figure l8. | 

j This means that for a certain jsuction amount, the influence range 
of a suction- slit becomes a maximum. For large suction amount coef- 

I ^ ' 

■ficients, the laminar boundary layer behind the slit becomes thinner, 
but not more stable. 

, Figure 19 shows the example of a Iplanej boundary layer. The shape 

f \ 

parameter inf-luence behind the suction slit decays faster, the thinner 

the boundary layer.-.ahead-.of the. slJ^t-.— - Thereof or-e-^- -we- can already draw 

'the conclusion that as the boundary- layer thickness increases, the slit 

^separations increase in the flow dlrectioh. 

: I 

‘ - I 

' i 

I 

'7 . 3 -‘3 Comparison Calculations 

j As in a test, example for the approximate solution for determining' 

■the laminar boundary layer influence by suction.-'.slits we can use a 
'suction strip over a plane plate, which has been discussed by several 
authors. | - 

1 . Krause [64] calculated the influencing of the boundary layer by 
suction strips of various widths for a constant amount of sucked-off 
‘flow, as shown in Figure 20. If inj the case of slit suction, the cor- 
responding suction .amount coefficient is 

^ i ^ , 

i -i •» _ s ~fRe I ' _ s 

• j ^ ''w Rej-2 0,6641 ' 

t 

. f 

; If the boundary layer influence plotted against the width of the /50 
suction strip is extrapolated to thL width of the suction slit, then 
we find the following: | 

- The thickness parameter reduction through a suction slit 

agrees well with the thickness parameter reduction through a suction 
btrip . 

, - The shape parameter increase (H-jo - ) through a suction 

slit can not be completely verified by the exact solution for a suction 

: 5 '8 

'strips for example, we have to consider the fact, that afTRe^^ = 10 , the-, 
narrowest investigated suction strip is about 220 times wider than the 
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corresponding suction slit. | 

The comparison calculation gii|en in Figure 21 shows that for the 
same suction amount, the influence on the boundary layer behind a suc- 
tion slit according to the amputation method and behind a suction strip 
according to the calculations of ijheinboldt [56]; Smith and Clutter 
[57]j Bethel [59], and Krause [64] are very .similar if the slit is ar- 
ranged in' the middle of the strip]. The relationship . 

.... _ 

, s/x V 

* o w 

! c = 

^ ° ^l+s/2x^' 0/6641 ,, 

exists. . -“71 771 ' 1 "-T"" ITT": . 

1 

The reliability of the amputat'ion method has therefore been proven. 
Therefore, the results of the slit suction calculation apply approxi- 
mately for narrow suction strips. | 


t 


7 . 4 Sink Effect' of Ring-Shaped Suction Slits _ 

An additional velocity is induced through each suction slit, which 
is superimposed on the potential flow. /51 



In the case:-of a body of revolution, we have the problem of esti- 
mating the additional velocities induced by ring--shaped .suction sinks'. 

j 1 

¥e will assume that a source ring Cxq, r^) with yield ► 
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= const.., 'induces the followin g a xial and radia l velocit y components.; 
a.t aT'poinir v') ' 


2Trr^ 


(7.24) 


_ 


Rr 


2^r, “qr 


The reduced, induced velocity components: 
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were tabulated by Kuechemann and .ifeber [ 65 ] as well as Dreger [66]. 

Ilf we introduce for the yield q for the suction amount coefficient cJ 

O CJ 

■ 3.C cording to(7.b)s. 


Q 


2ir I 


= _ 


(7.25) 


2irr„ 


then it is possible to write the Velocity components in a dimension- 
less form 



u 

CD 


= - 0/1^ 

(r„ /l)^ 


Uq 



u 

CO 


= . 0/1^ 

4ir^ (ryo^ 



(7.26) 


The tangential component— 




00 


'x dr T 


Qo 


dx u 


SSL. 




CO 


(7.27) 


is the additional velocity induced on the body surface by a suction 
ring. 

For small suction amount coefficients, the sink effect is restric- 
ted to the immediate vicinity of the suction sink, so that when there 
are several suction sinks, only the Influence of the two adjacent sinks 
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( 7 . 28 ) 
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CO 




u 

oo 


+ 



u 

oo 


has to he considered in the range 
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The additiona l v elocit y that is po sitive ahe ad of t he suction _ 
sink is negati ve b ehind th e sink . i means that there is a pres- 

sure drop due to the sink effect in the vicinity of the suction sink, 
which is superimposed on the shape-determined pressure distribution. 


= 1 -( 


ges 




TT.T97 


oo 


i i_ 

The sink influence is estimated for a circular cylinder in axial - 
flow having a suction ring in the center, which is given in Figure 
22. Because of r = r^ = const., and dr/dx = 0, we have the followingj 
for the circular cylinder 

1 

u u i 

_ oo CO 

Since the .slit: width in relationship to the boundary layer s/6 = 
1,1 the additional velocity in the immediate vicinity of the sue— 

-V 

tion slit depends substantailly on the boundary layer thickness and 

f 

Increases with increasing Incident Reynolds number. 

As can be seen from Figure 22, the additional velocities induced 
by ring-shaped slits over a body of revolution in the case of laminar 
boundary layer suction are so small compared to the Incident velocity 
over the Reynolds number range of interest that they can be ignored. 


7 . 5 Optimum Suction 


7.5-1 Optimum Position of Suction Slits 


In technical applications we are interested in the configuration 

of the suction slits, so that a body surface can be made laminar using 

the minimum number of suction point s.* /5^ 

^Technically we are also interested in the question of how the body surface 
can be made laminar with the minimum amount of sucked— in flow for a speci- 
fied position of the suction slits. The maximum possible slit separations 
are then taken from the present paper. 
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The suction relationships in Chapter 7-3 say nothing about the po- 
sition of the suction slits, but instead contain the following free 

( 

I parameters : 


- the boundary layer parameter ahead of the suction slits ;!Zj. , H 

- and the relative suction heights': n 

i iQi 


1 ■ 32i. 


; A condition for avoiding boundary layer transition is that the la- 

minar boundary layer remains stable i-n the entire suction:, region. If 
‘the state is partially unstable, the boundary layer perturbations would 
rbe buijlt] up in an uncontrollable manner in spite of suction, and this 

• would lead to premature transition! The stability criterion (-6.1) 
therefore is a boundary.—condl.t.ion .for optimum .suction of the .laminar 

'boundary layer. 

* The first suction slit is an exception to this, ahead of which one 
jean allow an unstable laminar boundary layer-, as experiments bv Pfenni- 
;nger have shown [6] - if there is a sufficiently long suction region 
•behind it, which serves as a damping path. Therefore, we-will introduce 
jthe transition criterion C6.8) as the boundary condition for the first 
(Suction slit. With this assimiption we can bring about a substantial 

' I ■ 

1 reduction in the number o-f slits, because as will be shown later on., 

I the slit separation required for the stability criterion is small. 

I By introducing the stability c^ondition or the transition condition 

(into the suction relationships, we [^ecify th^ boundary layer parla-par'; - 
meters ahead of the suction slits, f The position' of the suction slits 

» j ^ 

[depends only on -the relative suction height for a specified body 

t 

■shape and incident Reynolds number: 

t \ ^ - 

The region of influence of a suction slit can not be arbitrarily /55 
(extended by increasing the relative suction height, as we showed in 
Chapter 7-3.2. It becomes optimum !at a relative optimum suction height 


n 


Qopt ' 


When optimizing the relative suction height for maximum slit sepa- 
ration, we will assume the same relative suctionheight for ' all' s-uctlon 
slits . 


The end of the suction region 


is one of the free variables. 


*7-5.2 Total Sucked-In Flow 

The total suction amount of a jbody in a flow with n^ suction points 


la 
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The total suction amount coefficient of a body of revolution is 
calculated from 

c = Q = _^ I :l c 

° 0/1^ *=' ’ ' (7.30) 

where the local suction amount coefficient is defined by (7.6). 

Accordingly, the total suction amount coefficient of a plane con- 
tour is 


c 

Q 


Q 





(7.31) 


7 . 5.3 Consideration of Suction Power in the' Drag Balance 

When evaluating a body in a flow with boundary layer suction, it 
is necessary to consider the required suction power in the drag balance 
calculation. The references Edwards [67] and Torenbeek [68] contain 
approximate calculation methods for the suction power for special suc- 
tion installations over wings. 

Assuming that the suction amount does not contribute to the thrust, 
that is", it emerges with u = u^ and p = p<„^ at one point of the body 
again, the following pumping power is required for the suction slit 
>(i.) in the case of incompressible flow /56 


n N . = Q. 2. u - Q. (Ap. +Ap, s 

Qi I 2 00 I ^1 ^v.) 


C7.32) 


The pumping efficiency Pp and the pressure losses Ap^. on the suc- 
tion points and throttling points,, as well as in the lines, are only 
known after design of the suction installation has taken place. The 
publications of Pfenninger [69] and Gregory [63] state that the total 
pressure losses of a suction installation are rarely greater than 
APy = CQ.li+.(0.2) q^. In order to estimate the suction power, it is 
therefore possible to ignore the pressure losses, and we can set the 
pumping efficiency = 1 so that approximately we have the following 
for the suction point !i(l)(. 
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Q. (q -Ap.) (7.33) 

and for suction points, 

No = I Q| (q,^ - Ap.), (7. 34) 

Instead of the suction power, we will define an internal or equi- 
valent suction drag 



Wo = — 


(7.35) 


lAcq'ordingly , using C7.6) we can calculate the equivalent suction drag 
coefficient of a body of revolution 


W. 


2ir n 


q^O w S T Si 


C7.36) 


and of a plane contour 


Wo 




V == ^ 0“^ ) . 

Wq F Qi Pi 


(7-37) 


We then obtain the total drag of a body with consideration of the 
suction power from the sum of the external and internal drag values /57 


W = W, + , 


C7.38) 


Accordingly, the total drag coefficient of a body of revolution, 
with respect to the surface is 


W 


Wo 


q O 
^oo 


= c + c 


C7.39) 


and with respect to the volume, it is 

W 


O 


W, 


loo 


which is determined from C5.38] as c„ . 

In the plane case, the total drag coefficient is 


C7.^0) 
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W 

; = =- = c + c I 

Wc q r Wv, Wn . 

F ^OO “ I 


(T.itl) 


where is given by C5.39). ; 


7-5.4. Comparison Calculations 

As'a test example for a body in a flow with slit suction^ we can 
use the suction experiments of the Pfenninger group [6] with the test 
aircraft P-94A. In the test experiments the boundary layer was held 
laminar along the top side of a wing without sweep back using 69 suc- 
tion slits in the range 0.410 < xU /I < 0.953 up to the Reynolds nura- 

n ilD S 

ber Re^ = 3.6 x 10 . It was kept | laminar into the region of the trailing 
edge. j 


This lift coeffi- 
3 X lo"^, 


Figure 23 shows the contour and the pressure distributioncof the 
profile topside for the lift coefficient c^ = 0.15. 

cient corresponds to- the incident | Reynolds number Re^ = 3 x 10' , where 

the comparison calculations were made. 

The position. of the first suction slit in the experiments agrees 

very well with the calculated transition point positions for the case 

without suction, x., /I = 0.410. ' 

^ Abs . I 

Using the computer program discussed in Chapter 8j, however, we can 

not simultaneously maintain the suction conditions of Pfenninger 


- slit number n„ = 69 ’ 

^ 1 

- total suction amount coefficient c^ = 3 x 10 


_i} 


|b^*cause they are not optimum withi respect to minimum slit number. 

There fore in the comparison calculations we vary the parameter ^ so 

experiments each was satisfied. The 


'in, Figures 2 3 |andM2|4^ 


that one suction condition of the 
comparison calculations are given 

I 

- if n„ = 69 is specified, then the calculation gives about the 

y. I 

same slit separations, a substantially reduced suction amount, and a 
•somewhat reduced total drag compared with the' suction experiment.. 

- if c^ = 3 X 10“^ is specified, then the calculation gives sub- 
stantially less -suction opcints but about the same total drag. 

We should consider that the suction experiments were performed at 
the Mach number M = 0.6, but thatjthe calculation was restricted to ^ 
incompres-s-i-b-l-e— flew-.--- - ' - • - ‘ . 
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The comparison calculatiors confirm that- the equivalent suction drag 
I can be approximated using the appr|Oxlmations agreed to in Chapter 7- 5- 3* 

; Therefore 3 the drag coefficients calculated for slit suction are valid. 

' If the relative suction height is suitably selected, it is possible to 
: substantially reduce the number of slits and the suction amount compared 
with the suction'- 'configuration of Figure 9- 

/59 

1 8. Description' of Computer' Program 


8.1 Preliminary Eejiarks, . | . 

: A computer program was developed for characterizing the boundary 

* 

j layer and drag of' bodies in flows i|with laminar influencing suction 
1 slits. It allows the following sy^stematic variations: 

; - rotationally'-symmetric or plane cont.ours. 

’ 5f-thickness ratio d/1 1 

t 

i f Nose,i lg/1 


^length ratio of '{Middle part Im/^ 


*body function of) .. 


Tail W1 
i 

nose L = 


tail = fCCti) 


- the suction conditions for slit suction 

I 

* relative suction height 


” suction -region /I 


Beginning 


lEnd 


number of suction points] n 


- transition point position x^/1 in the case without suction 


* natural | 

completely laminar 

* forced, especially . / 

completely turbulent 


.- 3.2 
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- and the incident Reynolds number 
Using the program, the following are calculated: 

- the boundary- layer variables at the support points Z/1, ^-^2^ 

Re 2 '= f(Xg/l). 

-The boundary layer and suction variables at the suction points 

'*1.1 !'^'l ^321 - V‘- ^3 2, 

-The boundary layer variables at the points 

Z/l, H 22 , ^6^2 = f(x,^/l) . 

-The total suction removal coefYicient“c“ 

-And the drag coefficients c , c , c , c 

Wq ' Wo ' Wv ■ 

8 . 2 Structure of Computer Program 
8.2.1 Numerical Methods 

Boundary Layer Calculation Using the Step Formulas (5>22). 

The boundary layer variables Z, H ^2 ^.t the end of interval Ci) are 
determined by iteration from the values at the beginning of the inter- 
val (1-1). For the first interval iteration step v = 0, the variables 
at the previous interval (i-1) are used, The iteration accuracy is 

I H3^2!^ - I - 


'The; Iteration converges rapidly 


U.T. 




- 1 




< 0,03 


( 8 . 2 ) 


and 


32j ^32j_i 


< 0,003 


(8.3) 


In the program, there is a step selection which is automatic, star- 
ting with the distance to the next, support point 


, X 

1 , 1-1 _ s 


X. 


T “ T' 


(8.4) 
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The step is cut. in half 


1 , 1-1 


1 

2 


Ax} 




i,i-l 


I 


(8.5) 


jwhen the convergence -conditions with the selected step are not satis- 
fied , or if the accuracy cannot be | reached in v = 10 iteration steps. 
A maximum of u = 80 step halving operations 


IS. 


allowed. Investigations 
of the time optimization of the boundary layer calculation were performed 
by Otte [70]. 

Determination of .the.- Except.iohaI‘ Po-in.ts. 'in^-.a- Bo.undary Layer Calculation 

i — — — > ^ 

The following are exceptional (points for the boundary layer varia- 
tion • (G) : 

j - suction slits 


- instability point 

- transition .point 

- laminar separation point 

- turbulent separation point j 

- end of the region where the l^otationally-symmetric boundary layer 


, calculation- is valid 


,for which the criteria .('K) are specified. The calculation of these 
jpoints leads '.to determining the zeijoes of the function 


G 


(t) - Fk<T> = ° i 


C8.6) 


with the ‘required accuracy 




X. 

I 

T 

n 


< 10 ! 


C8.7) 


Since the function (8.6) is not available in analytic form, the /62 
t ■ i 

zeroes are calculated numerically by boxing an. 


; 8 . 2 . 2 Program Structure 


The computer program consists 
- the main parts 

* Initial calculations 

* boundary layer calculations 


of the following; 
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t 

*Initial calculation 
- and the routines I 

*suction slit | 

*Mangler transformation. ’ 

In the "suction slit" routine,' the suction ivariables are \ 

\ 

shown in Figures 12 to 15. * 

In the "Mangier" transformatioh^routine , the Mangier constant [L/If 
and the average shape parameter , which are requird for the transi- 

tion point determination. Formulas' (6;17) and C 6 .I 8 ) are used. 


8 . 2.3 Special Conventions ■ ' 

' " ■ " ' ■ • 

The possibilities of the computer program can be found in Chapter 
'8.1. Additionally 5 the following conventions are agreed upon: 

- there is boundary layer transition at the laminar separation 
pointy ' 


- there is no drag calculation' for turbulent separation 

- the drag calculation from the boundary layer variables is done 

as follows: J 

*At the point = 1/15' for rotationally-symmetric contours, 

*At the last support point for plane contours. 

- the last boundary .layer transition occurs as follows: 

“at the boundary layer ^ 2 /^ ^ I/I 6 for- a rotationally-symmetric 

'contour. 


*at the next-to-last support point for a plane contour. /63 

- the suction- region ends prematurely, 

1 

*for a forced boundary layer transition 

*as-soon as I .5237 > Hoo ; > 1.6250 applies .ahead of the suc- 

3'=-Ii 

tion point. j 

If the laminar boundary layer is still unstable behind the first 
jsuctlon -IS lit, then the suction height is increased until a stable state 
is reached. ' 


I 

8. 3 ALGOL-Program | 

The computer program is -written in ICL-ALGOL for the ICL 1900 

I X 

computer Installation o'f the Berlin Technical University.. Compared 
-^The computer program is available" 'from the Institute For Aircraft Design 
for T.U., Berlin. 
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with ALGOL decLarations , there are pnly slight modifications. 

The program has a memory requirement of nominally 13 3 250 locations 
in this computer. For 100 support! points and 100 suction points, this 
is expanded to about 17,500 locations. In calculations without suction, 
the memory requirement is substantially reduced, if the program varia- 
tion is used without the suction slit "routine". 

The calculation time primarily depends on the number of slits. 
Without suction, it is about 30 seconds and increases to about 7 mini utes; 
for about 100 suction slits. J 

Figures 3, 7, 8 , 23, and 2H show test calculations with the pro- 

t 

gram. 

For systematic contour variations, the program can be coupled with 
a pressure distribution program, which contains a suitable shape sys- 
tematic program and a support point interpolation. If the input and 

outputs of the programs are coordinated, then the input of the support 

i 

point values is reduced to certain] parameters which characterize the 

body shape. A program for pressure distribution has been given by 

Oehler for bodies of revolution [9]. 

J 

< 

i 

9 . Results of the' Variation Calculation - 

The computer program discussed in Chapter 8 affords-the possibility 
of extensive parameter studies. We will select a few to give an idea 

I 

about technical applications. ! 

9 . 1 Plane Plate: Variation of Suction Conditions 

Already by using the example of a plane plate in longitudinal 'flow 
(pressure gradient dp/dx = 0), we can obtain basic information about 
'the optimum configuration of the suction slits. In this case, the 
most Important inflluences] in the variables are the following; 

- the incident Reynolds number Re-, 

! j 

- the relative suction height ^tiq 

- and the suction region 

The Reynolds number Re'-^ is given by the -problem formulation. The 
' suction jhel^t tIq should be~| optimizedj for a minimum number of slitfe^ if 
we agree that the suction region is such that .it starts in the transition 
point and ends for a relative plate length of x/1 = 0..95- 
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Figure 25 shows the boundary jlayer variation for Re^ = 10 and 

ri^ = 0,1. Figure 26 shows the suction variables plotted over the 
‘ 1 

plate length. As already discussed in Chapters 7*3.2, the suction slits 

with small suction heights bring allout a substantial shape parameter 

increase and a slight thickness parameter decrease. The boundary layer /65 

‘thickness increases in the ■ flow direction, nevertheless. However, the 

suction action behind a slit: decreases more slowly the thicker the boun- 

'dary layer is. At n„' = const, the slit separations become greater with 
' y 3 _ i 

'path length. Accordingly, the optimum slit widths and the local suction 

^amount coefficients Increase in thi flow direction. In the literature,, 

' ) 

this basic fact has,- not. yet., been. p-uhli shed... - 

, I 

Figure 27 shows tha't the slit {separations, slit widths, and the 
local suction amount coefficients are reduced within -increasing inci- 
,dent Reynolds number because of the decrease in the boundary layer 
.thickness. • i 

’ The increase in the number of Jslits with.' increasing inc-ident Rey- 

jnolds number shown in Figure 28 is |caused by the earliest beginning of 

^suction and the reduction of the sljit separations.. The reduction in 

|drag due -to suction is also considerable, if one considers the suction 

'power. For = 0.1, the drag coefficient, of the plate which has been 

I y ' I 

-made laminar by the suction slits is only slightly above the Blasius ' 
curve. The total suction amount coefficient first increases because 
iof the increase in the number of 'sl>its with incident Reynolds number. 

It then reaches a maximum value, and then decreases again at high Rey- 

^ i 

nolds numbers, because then the red^uctlon in the suction. |^o_unt coeffi- _< j 
{cients has a greater effect than the Increase in the number of slits. 

If the suction already starts at the Instability point, then the slit 
number increases at all .Reynolds numbers. The drag coefficients and 
^the suction amount coef flcients '.onl^ increase in the- lower Reynolds 
number range. 

i ‘ - ^ 

, The results of the simultaneous variation of relative suction height 
and infinite Reynolds number are sh^own in Figures {29 and 30. The ^ 
optimum .relative suction height is almost independent of the Reynolds 
number, at which the number of slitjs becomes a minimum, 

' ' i Iq opt Iq min ^ ; 

If less air is sucked away at jthese suction points, then ihe drag ’ 
and. the.. .totah_.s.uciipn^amQnnt_d^jcrease.,,„hu.t...th,e,=xe 9 .uired' slit mimber 
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increases rapidly (tIq < 0.1). At = 0.02, one comes very close to 
the case of continuous suction, see Figure 29. An increased suction Pq> 0.1 
is out of the questioinji then the slit number and the drag, as well as the 
total suction amount, increase considerably. 


9 . 2 Bodies of Revolution 

9.2.1 Variation of Suction Conditions 

For a body of revolution we will investigate the ways in which the 
influence of the pressure gradient iaff ects the optlmdm c onfiguration ” 
of the suction slits, compared with the plane plate. For this purpose, 
we must vary the incident Reynolds number Re^, the relative suction height 
tIq, as well as the suction range x^^^/1. 

As far as selecting the body shape is concerned, it should be men- 
tioned that without suction, spindle shapes with a pointed nose part have 
been suggested, see Hertel [1 to 3]. These satisfy the condition for 
minimum specific drag much better than conventional shapes with a long 
cylindrical central part and a small thickness ratio. It has been found 
that these spindle shapes are favorable in the -case" of suction, compared 
with conventional shapes.' 

Since the nose contour of low drag spindle bodies is approximately 
parabolic according to [2] and [3]j as a basic shape we will select the 
following symmetrical paraboloid for the variation calculation, where 


the length ratios are 



the thickness ratio is , I ' „ 

1b = Sb 

This basic shape can be matched to various problems by changing the 
thickness setback and the thickness, ratios ^ From the transition and /67 
drag curve of the basic shape shown in Figure 31, we can see that the 
"laminar effect" is lost at high Reynolds numbers. As we will now 
show, making the flow laminar using slit suction in this area leads to 
substantial drag savings. , 

Figures 32 and 33 show the boundary layer variation and the suc- 
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it Ion variahlesi- of the paraboloid distributed over the body shape, for 

{ 7 

the incident Reynolds number Re^ = 5 x 10 and the suction conditions 
’Pq = 0.1 and = 0.5- They are very similar to those of the plane 

: plate shoxNrn in Figures 25 and 26. f The length of the laminar buildup path 
: behind the last suction slit is substantial. 

Figure 3^ shows that in-spite|of the reduction. in the stabilizing 
■pressure gradient, ja^ p„ = const., {the slit separations of a spindle 

^ L_iJ y I 

body are increased in the flow direction with the exception of the tail 
'end, where there is a steep pressure increase.- The relative slit width 
' and the local suction amount coefficient depend on the Reynolds number 
and Increase greatly .because ^of the -thickening- -of -the boundary layer at 
;the tall end. -Increasing the slit f separations in the tail area of the 
[bodies of revolution comes about because in the case of bodies of revo- 

j 

jlution, the reducing pressure gradient is opposed by an increased boun- 
; dary layer thickness, due to a decrease in the body radius. This effect 
fdoes not occur in plane flow, as shown in Figure 23. . - 

I For applications, various suction regions are possible for bodies 

jof revolution, for example, • . 

[ - in the case of ' an aircraft body, -it only seems that keeping the 

[nose part of the aircraft laminar is promising because of the disturbing 

I ^ • 

1 influence of the 'wings, . . | 

i i 

) - it seems promising to make the entire body surface of an under— 

' i 

(Water vehicle laminar using suction slits. 

1 Figures ■ 35 and 3^ show the results of a suction region variation 

; and Reynolds number variation at nJ = 0,,1. ■ Already a small suction 
I - y 

;region with a corresponding low suction amount is sufficient to -keep 

the nose part laminar, and this is jassoclated with a substantial savings 

,in drag. The entire body surface can be kept laminar with a small - /68 

'number of additional suction slits , but the total amount of sucked air 

[increases drastically. The extremely low drag coefficients in the <case 

iof suction to the area of. the tail tip are caused by the fact that when 

;the body surface is made completely laminar not only is the friction drag 

'-reduced--, but the pressure drag is also greatly reduced. The jump in the 
I j 

.displacement of the transition point, and the reductioniin the drag coef- 
ficient due to a suction slit are covered by special conventions in the 
{computer program, see chapter 8,2.3. 

' When there is a simultaneous variation of the relative suction 


59 



height and. the incident Reynolds number. Figures 37 and 38 show that 
the- slit number becomes a minimum If or bodies of revolution, indepen- 


dent of the Reynolds number. 

; 'lo 


If -t|he suclion height is 


'Iq 


opt 


0,1 


Figure 39 shows .that the optimum relative suction height is also approx- 
imately independent of the suction ‘region: The" variation of the suc- 

tion conditions over a .spindle body shows that = 0.1 is an optimum 

suction parameter for technical applications of slit suction. 

5 

I 


i 


9 . 2.2 Variation, i:Of Body' Shape 


In applications we are interested in the effect of the change in 


the body shape of a body of revolution on' the configuration of the 
suction slits. 


' At hQQp.j- ="0.1, we will investigate the influence of 

i - the thickness |setback .j or th|e relative noselength (5 b/15 

I - and the thickness ratio d/1 ' 

1 f 

I on the minimum .slit number. ■ We will start with the base shape dis- 
cussed in the previous chapter. | 

f Figure 1|0 shows the contour and' the pressure distribution of spin- /69 

I I 

; die bodies where the thickness setback was varied with the same thick- 
! ness ratio d/.l = 0.2, and shows that in the case’ without suction, the 
[ bodies with a large thickness setback have the lowest drag in the low 


Reynolds number range, belcaus'e] of the long laminar buildup paths. At 
the higher Reynolds numbers, this tendency is reversed. The influence 
of the thickness setback on the optimum suction configuration CiIq = 0.1) 
is given in Figures .42 and 43." ¥j|th the exception of small Reynolds 
numbers, the slit number is 'substantially reduced with decreasing 
thickness setback. 

Figure 44 shows the contour and the pressure distribution of spin- 
dle bodies, in which the thickness, ratio has been varied for the same 
thickness setback P]g/p| = Q.5- Figjure 45 shows that in the case without 
suction, the thicker bodies have a^ higher drag. One exception is the 
medium Reynolds number range, wher^e the long or laminar buildup paths 
produce a smaller drag for the thijcker shapes. It is remarkable that 
for a thickness ratio d/1 = O. 3 , there is still no turbulent boundary 
layer -separa.tion'..at_the._t.all... _ ■' _ 
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: The effect- of the thickness yariation on the optimum suction -con- 

figuration .(n^ = 0.1). is shown in Figures 46 and 47 . As the body thick- 

-w I 

• ness increases, the slit number- decreases in the entire Reynolds number 
range . 

• The thickness setback variation and the thickness variation of a 


' spindle body clearly show that the 


pressure gradient over the nose part 


.is most important for determining the minimum number of- suction points 
.even if suction takes place near the tail tip. 


10. Summary . - . . 

: t 

: In aviation and underwater theory, only slits, gaps, or perforated 

i strips^ can be used for practical boundary layer suction applications. /70 
.Up to now, continuous suction has been almost the exclusive case which 
;has been treated, since the investigations are usually 'restricted -to 
plane contours . 

I In the present paper, we discuss the technically interesting case' 

j of keeping a boundary layer laminar using suction slits over bodies of 
[revolution in incompressible flow. |_ 

I The calculation of the 'rotatimally-symmetric laminar and turbulent 

‘■bo,undary layers is done by the integral conditions for momentum and 
'energy, and assuming a single parameter velocity profile. 

! A new serai -empirical criterion “is introduced for determining the 

i laminar-turbulent transition- point jof bodies of revolution, and its 


reliability has been demonstrated from transition measurements. 

1 

The changed boundary layer variables behind' the suction slit are 
[determined approximately' from the remaining profile which remains after 
.suction. By comparison with exact [solutions-, which are available .in 
the special case for a. narrow suction strip, we confirm the reliability 
•of this method. Prom these results we draw the conclusion that the 
, effects which can be brought about 


by the (suction j -slits can als o be 


approximately- reach ed by perforated suction strips. 


configuration, we calculated, the 
This calculation was 


In- order to evaluate the slit 
drag with consideration of the suction power. 

'confirmed by suction experiments. 

The ALGOL program was developed for the systematic variation of 
body shape. Incident Reynolds number, and suction conditions which also 
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includes the case of no suction and the calculation of plane contours. 


In addition to the stability ,and transition conditions, we also 
iintroduce the relative suction height and the suction region as para- 
'meters in the suction relationships. For constant relative suction 
(height, the slit separations in the flow directions are increased, as 
'•long as there is no great pressure increase. 

' The local suction amount coefficients and 'relative slit widths 

\ 

.'increase greatly .towards the tail tip. 

( ] 

The variational calculations show that, independent of the body 


1 shape, incident Reynolds number and suction region for the relative 


^suction height, th.e .slit- separations .are. .a.,maximum. and consequently 
,the number of slits are a minimum. [ Compared with the continuous sue- 

. i 

[tion Cnn ■^ 0 } 5 tihe suction [amount coefficient increases dras— 

:tically at = 0.1. 'The drag coei^ficient only increases slightly 
[with consideration. of suction power. 

j Finally, we Investigated the influence of the body of revolution 

shape on the optimum configurationjdf 'the suction slits. It is found 
that pointed spindle shapes with, a jlarge thickness ratio and a small 
(thickness setback can be made laminar using a relatively small number 
lof suction slits due to the large pressure drop over the nose, even at 
I high Reynolds numbers. | 

I The combination of favorable bpdy shapes and optimum suction con- 

jditions jleadsj to solutions which satisfy the condition for a small 
(flow resistance with a small number of suction points and a small suc- 


,tion amount, for the most part, at high Reynolds numbers . 

; . i 
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Figures 


Body shape r(x) 



7 





approximation: according to Eppler [23] 

- - - according to Walz [11] 
according to Geropp [l6] 


^plane 

, ' stagnation! 



Figure 2: Relationships for laminar velocity profiles with suction 

influence . 
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Comparison; - Delayed stagnation point flow 
Laminar separation point 

------ Exact solution of Howarth [24] = 

Present calculation method with approximations of Eppler 

1^23] x*,yi-Cl1196 



Figure 3 : Delayed stagnation point flow. Laminar separation point. 

Calculation method - exact solution 
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stability Calculations 
without suction 


w'ith suction 


, ° Tollmien, Lin (Blasius Profile) ; ^;Bussmann, Muenz (asymptotic suction 
[ 31 ], [ 32 ] profile) [ 36 ] 

Schlichtingj Ulrich (Pohlhausen i_i Ulrich (Iglisch Profile) [37] 

P6 Profile) [33] 

_^;Pretsch (Hartree Profile) [3^] l.* Ulrich (Schlichting-Bussmann 

profile) [ 38 ] 

Tetervin (Hartree Profile) [ 35 ] 


Stability Criterion : 

Approximation i«s ot,, . «,556 -76.«7n.67o-H„)'-’‘' 



Figure Stability of the laminar boundary layer with pressure gra- 

dients and suction. 
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ty profiles according 


Transition Measurement 


Wind Tunnel 


Free Plight 


• Jon*a [48] 
a Zalovctk, Shook [ 49 ] 

A Setaore, Zalovcik, Platt [50] 
T Oaviea, Saith, Hfgton [ 51 ] [52] 
^ Plascott, Kigton, 

Salth, Braoaell [ 55 ] 


Rotationally a SoltZi Kenyon, Allen [55] 

Symmetric - ■ 

Plow 


Plane o Brasloa, Visconti [ 47 ] 
Plow □ V. Oocnhoff [46] 

A Schubauer, Skramstad [30], 
V Boltz, Kenyon, Allen [ 54 ] 


Transition Criteria 
^Approximation 

1,515,* Rjj =* 1,560 : log (Retj^ -Re^jj ) = 1,6435 - 24, 20 ( 1,51 50 - Rjj ) 
1,560< H 3 J 1,625 : i og - Retjj ) = 3^51 2 - 967,5 ( 1,6250 - Hjj 

Buildup a = e“^ 

_ _ _ _ Empirical 3 according to Granville E^2] 


lllartree separation *Blasius Profile ■ j plane stagnation 



Figure 6 : Length of the buildup paths for natural transition of the 

laminar boundary layer with pressure gradient. 




Comparison: Plat Plate 

k 

- - - Known solutions 
Present computation method 



Figure 7: Plat plate in longitudinal flow. Instability and transition 

point positions.' Drag coefficient. Computation method - 
known solutions. 
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Comparison: 


o }’ measurement of Boltz, Kenyon, Allen [5^, 553 

present calculation method 

- - - present calculation method 


0^1 - 
<D 

^ d/i. auo — -^2 
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1 measurement / 

V 1 / 

4 ) ' 

" ^ 

in 


^ -a.2- 





d/< . a 1 3J - 


measurement 


A 


1»0 


_ P CV2_ OiQ 1i0_0 dj <;< d6 08 

^lative iDody lengthx/£l 1 re lat ive "body lengthix/& 



Figure 8: Body of revolution. Transition point position. 

Coefficient. Calculation Method - Measurements. 


Drag 
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’suction slit 





suction strip' 



; suction gap' 


Figure 9; Possibilities of local boundary layer suction 


— f 



Figure 10: Ratio of Boundary Layer Thicknesses^ 


Geropp [l6J model. 




j relative wall separation! r,®= y/&' 


Parameter: Relative suction height \ 



Figure 11: Laminar Velocity Profile ahead of/behlnd a suction slit. 

Single parameter model of Geropp [1^]. 
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Figures 12-15: Influencing of laminar boundary layer through a suction 
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Figure l4; Nondimensional Suction Amount Coefficient 


Figure 15: 




CO 
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Ahead of 

suction slit: Blasius profile ’.s’zs 




Figure l6; c 
Figure 16-17 
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Mna)' 


Figure 17; |H32^ , Z^/Zj , s* = f (cj) 
Influencing of the laminar plate boundary layer by means of a suction slit. 
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Parameters: nondimensional suction amount 



I ~ - ' nondimensional 

'suction amount coefficient -at 

I I Parameters: thickness parameter ahead 

p I of suction slit compared 

(Z!„ /x„-a« 09 ) I vith the value for a flat 
suction slit i plate 
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Figure l8: 


Influencing of Boundary Layer as a 
Function of Suction Amount Coeffi- 
cient . 


Figure 19 ; 


Influencing of Boundary Layer 
as a Function of Boundary Layer 
Thickness ahead of slit. 


Figures 18-19: Influencing of laminar plate boundary layer by means of a suction slit. 
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of sucked-ln flow 

Figures 20-21; Influencing of laminar plate boundary layer through ^a suction point. 
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,1Modified Profile Topside^ 
MCA 65 - 213 I 





Figure 23: Position of suction slits. Figure 24: Equivalent suction resistance coef- 

ficient c^Q. Drag coefficient 

. r, T 1 4-4 with respect to the suction power c,™ 

Figures 23—24; Wind profile with slit suction. Calculation ^ wP 

method - measurements of Pfenninger [ 6 ]. 
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figure 25: Boundary Layer Variation 



Figure 26; Local Suction Variables 
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Figures 25-26: Flat Plate with Slit Suction 
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height] parameters : Incident i 

Reynolds number jxr, 



Figure 27; Slit separation, slit width. 

Local suction amount coefficient 

Figures 27-28; Flat Plate with Slit Suction 



Figure 28; Slit number, drag coefficient, 

total suction amount coefficient 
for various suction beginnings. 


REPEODUCIBirJTy OF TH^ 
0BIOL-.AL PAGE IS POOR_ 






Parameters : 

Relative suction height 


A 

OJ 

rS 


End of Suction Region 


I Parameter : 

I Incident Reynolds numher Re^ 



o ^<1 


, , „ without' 

!p suction^ 

§ ro o cvoooe-^s^transition 

•H U - - - --. 

•H C! I CV COOS'! 

tH *H O ' 

(i0004i 
O ' 

j3(Ci000Z 

m , 

0 

^ -d ,Ql 


^05 suction^ 


<H W 

(u a 
o o 
o o 





Tncideht Ueyholds number] r«, 


100 


e . 

, _iO o. 
•P -H ‘ 

^ -H 
o cd 

fH 
0 ) 

•H 

W 

« 

O 
O 


C^ooio 

Is 
o 
ft 


bo ^ 
cd P> 
(H 'rj 
W > 


Q0006-I 

C^0004 

dOOOZ 

0 


Tend of 

Tsuction range- 
»*ME A -0,9$ 







. 

___J_ 10 h— 

— 

10 * — 





1405 CV10 0H5 
[relative suctit)F’hieigEt“ i 


0^20 


Figure 29 : n 


Q" wO" 


= fCReO, Parameter ri 


Q 


Figure 30: n 
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Figure 29-30: 


Plat plate with Slit Suction, 
coefficient ^ 


Slit number, drag coefficient, suction amount 
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Figure 33: Local Suction Variables 
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Figures 32-33: Bodies of Revolution with Slit Suction. 




Slit Width. Local Suction Amount Coefficient. 
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Parameters : I 

end of suction I 

• I region 

Relative suction height'! 



Parameter / ^ 


[relative 
\ suction height" 

io-ai 


Parameter; _ 

I Inc'ident~Re^olds Number «», 






Parameter Re^ 


Figures 35-36: Bodies of Revolution with Slit Suction. Slit Number, Transition Point Position, 

Drag Coefficient, Suction Amount Coefficient. 
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Pigures 37-39: Bodies of Revolution with Slit Suction. Slit number, drag coefficient, suction 

amount coefficient. 


8 


REPEODUCIBILITY OF THE 
OEIGMAL PAGE IS POOR 








Parameter : 

Eelative nose lengths 


r- surface ratid"' 

d// = Q 2 



(a/t d/J = Q2 



'Parameter': 

Relative nose lengths' 



Figure 40: Contour and Pressure Distribution Figure 4l; Transition point position and drag 

coefficient. 
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Figures 40-4l; Bodies of revolution witbi various thickness setbacks. 
















Par a mete r: Incident Reynolds numter . 





relative nose lengths-, Js/i 



relative nose lengths , 


ELgure 43: 


Bodies of Revolution of Various Thickness Setback.s with 
Slit Suction. Thickness ratio d/1 = 0.2. 

no , cwo , Cq = f a i Parameter Re ^ • 
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Parameters tj - nose length .|/b/' = 5 
thickness ratio.] 0, c "tail” length] <h- ■f = as 


^Parameters \ nose length (g/.- - 45 
jthickness fati^' <s/i [tail length <h/< - as 
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Figure 44: Contour and pressure distribution Figure 45: Transition point position and drag 

point coefficient. 


Figures 44-45: Bodies of Revolution of Various Thicknesses 















['' Paramet e r : In cident Reynolds number. Re ^ 

|_enA suction region] '^end oT"' sucti orT on ; 

=‘Ab.E/^=^5 ><Ab.E/^=Q3 



Figure 47: Bodies of Revolution of Various Thicknesses with Slit 

Suction. Relative nose length lg/1 = 0.5- 

Hq , Cv,o , Cq = f ( d/i). Parameter Re^ , ^ • 
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: 13. Appendix 

: Calculation of Integral Expressions use'd In Chapter 7-3»l 
’ The integral expressions ( 7 .I 8 ) and ( 7 • 2 3 ) 
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are calcu.l^ed by introducing the following coefficients 


Csee C7*12) 
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•[(c + 4)(c + 3)(c + 2) + 

+ CL, (c + 4)(c + 3) + 

+ 2(32 (c + 4) + 
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The integral expressions in condensed notation are: 


_ doc) + (c + 4} 2 q 2 + (c + 3) a-| + (c + 2) 
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